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We compile from literature the broadband SEDs of twelve TeV blazars observed simulta- 
neously or quasi-simultaneously with Fermi/LAT and other instruments. Two SEDs are 
available for each of the objects and the state is identified as a low or high state according 
to its flux density at GeV/TeV band. The observed SEDs of BL Lac objects (BL Lacs) are 
fitted well with the synchrotron + synchrotron-self-Compton (syn+SSC) model, whereas 
the SEDs of the two flat spectrum radio quasars (FSRQs) need to include the contri- 
butions of external Compton scattering. In this scenario, it is found that the Dopplcr 
factor 5 of FSRQs is smaller than that of BL Lacs, but the magnetic field strength B of 
FSRQs is larger than that of BL Lacs. The increase of the peak frequency of the SEDs is 
accompanied with the increase of the flux for the individual sources, which seems oppo- 
site to the observational phenomena of the blazar sequence. We refer this phenomenon 
to blazar anti-sequence of spectral variability for individual TeV blazars. However, both 
the blazar sequence from FSRQs to BL Lacs and blazar anti-sequence of the spectral 
variability from low state to high state are accompanied by an increase of the break 
Lorentz factor of the electron's spectrum -yy, and a decrease of B. We propose a model 
in which the mass accretion rate M is the driving force behind both the blazar sequence 
for ensembles of blazars and the blazar anti-sequence for individual blazars. Specifically 
we suggest that the differences in (M) of different blazars produce the observed blazar 
sequence, but AM in each blazar results in the observed blazar anti-sequence. 

Keywords: Radiation mechanisms: non-thermal; BL Lacertae objects: general; quasars: 
general; gamma-rays: theory 



1. Introduction 

Blazars, a subsample of active galactic nuclei (AGNs), are composed of BL Lac 
objects (BL Lacs) and flat spectrum radio quasars (FSRQs). It is well known that 
the observed emission from blazars is jet dominated and the observed broadband 
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spectral energy distributions (SEDs) are bimodal. Generally, the bump at the IR- 
optical-UV band is explained with the synchrotron process ol relativistic electrons 
accelerated in the jets and the bump at the GeV-TeV gamma-ray band is due to 
the inverse Compton (IC) scattering ol the same electron population. The seed 
photon fiel ds m ay be from the synchrotron radiation themselves (the so-called 
SSC model or from external radiation fields (EC), such as the broad line re- 
gion (BLR)^. Broadband SEDs obtained simultaneously or quasi-simultaneously 
are critical to investigate the radiation mechanisms and the physical properties of 
the emitting regions for the two kinds of blazars. 

With a large sample of different types of blazars, Fossati et al. (1998) reported a 
spectral sequence ol FSRQ-LBL-HBLpI, i.e., along with this sequence, observation- 
ally, an increase in the peak frequency ol synchrotron radiation (i/ s ) corresponds to 
the decreases ol bolometric luminosity and the ratio (Lie/ L syn ) ol the luminosities 
lor the high- and low-energy spectral components. Theoretically, an anti-correlation 
between the break Lorentz factor of electrons (7b) and the energy density of radi- 
ation regions (U ) in the comoving frame was found. So this blazar sequence was 
interpreted by Ghisellini et al. (1998) with the cooling ol the external photon fields 
(such as the BLRs) of the blazars^. The discoveries of 'blue quasars' posed a chal- 
lenge to the previous interpretation for blazar sequence. To address this problem, 
Ghisellini & Tavecchio (2008) reported that, more physically, the sequence is due to 
the different black hole (BH) masses and accretion rates ol the sources^. According 
to this interpretation, the 'blue quasars' should have large BH masses and interme- 
diate accretion rates and thus the emission region is beyond the BLR. Moreover, the 
red low-luminosity blazars should exist, which have small BH masses and relatively 
large accretion rates. More recently, Chen & Bai (2011) extended this sequence to 
narrow line Seylert 1 galaxies^, which are similar to the low-peak-frequency and 
low-luminosity blazars. 

Observations show that the SED peak frequency ol a TeV blazar increases with 
its flux, indicating a tendency that a brighter TeV emission corresponds to a harder 
spectrum for the emission at the X-ray and gamma-ray bandd^M This phenomenon 
seems to be opposite to the blazar sequence; we thus refer this phenomenon to blazar 
anti-sequence of spectral variability for individual TeV blazars. The spectral shift 
at different states of TeV emission for the individual sources may be due to the 
different 7b at different states^. Thus, it seems that both the blazar sequence and 
anti-sequence are caused by the change ol 7b. However it remains unclear whether 
there are any connections between the two phenomena. So lar, about forty AGNs 
have been detected at the TeV gamma-ray band. Except for two radio galaxies, all 
of the confirmed TeV AGNs are blazars. Moreover, most of the confirmed sources by 
Fermi/LAT are also blazars. The abundant observation data provide an excellent 
opportunity to reveal the physical connections between the blazar sequence and 
anti-sequence. 
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2. Sample Selection 

In order to obtain well-sampled SEDs, only the TeV blazars, which have positive 
Fermi/LAT detections and two or more observed broadband SEDs, are considered. 
Twelve blazars are included in our sample, two FSRQs (3C 279 and PKS 1510- 
089) and ten BL Lacs (W Com, Mkn 421, Mkn 501, PKS 2155-304, 1ES 1101-232, 
BL Lacertae, 1ES 2344+514, 1ES 1959+650, PKS 2005-489, S5 0716+714). We 
compile their broadband SEDs that were simultaneously or quasi-simultaneously 
observed with Fermi/LAT and other instruments from literature. The two well- 
sampled SEDs for a source are identified as a low or high state according to the 
observed or extrapolated flux density at 1 TeV, except for PKS 1510-089. Because 
the SEDs of PKS 1510-089 are cut off at 1 TeV, the high and low states of this source 
are defined by the flux density at 10 GeV. The observ ation data of the broadband 
SEDs for 3C 279 and PKS 1510-089 are from literaturfPHHIl. 



3. Models and Results 

The observed SEDs for those sources are double-peaked. For BL Lacs the external 
photon fields from their BLRs, if exist, are very weak and negligible, compared 
with the synchrotron radiation photon fields. Therefore, only the syn+SSC model 
is considered to fit the observed broadband SEDs of the ten BL Lacs. The detailed 
calculation process for the model and the strategy for parameter constraints can 
be found in RefP Following the methods described in RefP, the SEDs for the ten 
BL Lacs are fitted well by the single-zone syn+SSC model; the fitting results and 
parameters can be found in Fig. 1 and table 1 of RefP. 

For the two FSRQs, the contributions of external field photons from their BLRs 
need to be considered. The total luminosities from the BLRs for the two FSRQs 
are taken from RefP and then the BLR sizes are calculated using the BLR lu- 
minosity with formula (23) given in RefP. The radiation from a BLR is assumed 
to be a blackbody spectrum and the corresponding energy densities seen in the 
comoving frame are L/blr — 5.37 x 10~ 3 r 2 erg cm -3 and Z7 BLR = 4.98 x 10~ 3 F 2 
erg cm -3 for 3C 279 and PKS 1510-089, respectively, where we take T = 6. Thus, 
the syn+SSC+EC/BLR model is used to fit the observed SEDs of the two FSRQs 
and the model can explain the observed SEDs well as shown in Fig. 1. 

The range of 5 for these BL Lacs is from 8 to 50 and is clustered at 11 ~ 30. 
On average, <5 of BL Lacs is larger than that of the FSRQs, but the magnetic field 
strength B of BL Lacs is smal ler than that of the FSRQs, consistent with the results 
of some previous work a^^ l The values of 7b for these BL Lacs vary from 10 3 to 
10 6 , but the values of B are clustered around 0.1 ~ 0.6 G. No correlation between 
7b and B for these BL Lacs is found. 
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4. The Blazar Anti-Sequence and the Blazar Sequence 

4.1. The Blazar Anti- Sequence 

As shown in Fig. 1 and Fig. 1 in Ref!^, there are clear spectral shifts accompanying 
the flux variations for these TeV blazars. The SEDs of the twelve TeV blazars in 
our sample are obtained in the low and high TeV states, which is defined according 
to the flux density at 1 TeV except for PKS 1510-089, for which the state is defined 
according to its flux density at 10 GeV. We compare v s and v c between the high 
and low states in Fig. 2(a). It is clear that the SEDs in the high states shift to a 
higher energy band than that in the low state. 

As reported in our previous papeJ^, to investigate what may be responsible for 
the spectral shift in the low and high states, we derive the ratios of the flux density 
at 1 TeV (i?i Tcv) and the physical parameters (R x ) in the high state to that in the 
low state for these sources, where x stands for Lboi, B, 5, 7b, or Pj G t- The results on 
the SEDs of PKS 1510-089 are not included here. The values of 7b, L^oh and -Pjet 
of the high states are systematically higher than that of the low state. A tentative 
correlation between i? 7b and R\ ToV is found with a correlation coefficient r = 0.55 
and a chance probability p — 0.077 as shown in Fig. 2(b). Therefore, it is possible 
that the spectral shift at different states is due to the different 7b of these sources. 

In order to compare the differences of parameters for the high and low states, we 
also calculate the magnetic field energy density {Ub), the available photon energy 
density (U ph , U ph = U BLR + U' syn ^ aa for FSRQs and U ph = U' syn ^ n for BL Lacs) 
for IC process in the comoving frame, the luminosities of the synchrotron radiation 
(isyn) and the IC process (Lie, Lie = -^ssc f° r BL Lacs and Lie — ^ssc + £ec for 
FSRQs) in high and low states. The ratios of those parameters for the two states 
(Rub-i Rtj' ' ^i.yni and Rl 1c ) as a function of the ratio of 7b (-R 7b ) are shown in 

ph 

Fig. 3. No correlations between the ratios of those parameters and the ratio of 7b 
are found. As presented in Fig. 3(a), for most of the sources Ub becomes smaller in 
the high state than in the low state. However, the ratios of U b , L syn , and Lie are 
larger than unity for most of the sources as presented in Fig. 3(a), (b). 

4.2. The Blazar Sequence 

As described in section 1, according to the observational phenomenon of the blazar 
sequence^, one can expect that both Lboi and Lie/L syn are anti-correlated with 
v s . Lboi and Lie/L syn as a function of v s , are shown in Fig. 4(a) and (b). A weak 
correlation is found for v s — £boi- The Spearman correlation analysis yields a cor- 
relation coefficient r = —0.58 and a chance probability p = 0.05 for the high state 
data, r = —0.61 and p = 0.04 for low state data, respectively. Excluding the two 
FSRQs, no correlation would be found. However, the ratio of Lie/L s is indeed anti- 
correlated with i> s , especially in their low states. The Spearman correlation analysis 
yields a correlation coefficient r = —0.71 and a chance probability p — 0.009 for 
the high state data, r — —0.94 and p < 10 -4 for the low state data, respectively. 
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Because most of the sources in our sample are BL Lacs, the external photon fields 
outside their jets are much weaker than the synchrotron radiation photon field and 
the EC process is thus not considered for these sources. As v s increases, the SEDs 
shift to the higher frequency end and the KN effect should be more significant. 
According to equation (20) in RefP, t/ s ' yn>avail = EC(f^) 1_Ql , the available 
photon energy density of synchrotron radiation for IC process decreases with the 

increase of 7b- For ^r 2 - ~ ^"t}^""" , the ratio of issc/^syn would decrease as v s 

increases, since Ub is almost constant and U syn avail decrease along with the increase 
of 7b for the BL Lacs in our sample. Therefore, the anti-correlation of L\q / L syn — v s 
may be also due to the KN effect, especially for BL Lacs. 

Because the interpretation for blazar sequence is cooling of the external photon 
fields, a more "theoretical" scenario than the purely phenomenological sequence is 
the anti-correlation between the break Lorentz factor o f elec trons 7b and the energy 
density of radiation regions U in the comoving frame^H^. 7b as a function of U 
V = ^ y „,avaii + U B for BL Lacs, if = U' syn ^ Q + U BW + U B for FSRQs) is 
also shown in Fig. 4(c). Although the correlation between 7b and U of our sample 
sources has large scatters, especially for BL Lacs, comparing our results with that 
m RefR the two results are consistent and the BL Lacs included in our sample 
distribute in the left top of Fig. 4(c), where the EC process is not important. If 
the different 7b is totally due to the different external photon field, there should 
be a correlation between the luminosity of the BLR and the peak frequency of 
synchrotron emission. With a FSRQs sample^l, however, no correlation between 
-t'BLR and u s is found as shown in Fig. 5(a). Sometimes, for simplicity, the luminosity 
and the radius of a BLR is assumed as Lblr = O.lidisk and -Rblr = 10 17 ^d(sk 45- 
So the comoving energy density of a BLR is given by Z7 BLR = 3.76 x 10 _2 r 2 erg 
cm~ 3 under these assumptions^! and is totally decided by the value of V = 5. 
Nevertheless, no correlation between T 2 and 7b is found either, as presented in 
Fig. 5(b). 7b is correlated with the total energy density of the emitting regions in 
the comoving frame as shown in Fig. 4(c), but not correlated with the energy density 
of the BLR. Using the sample data in RefP 

and our sample data, it is found that 
7b is correlated with the magnetic field energy density, as shown in Fig. 5(c). Some 
works indeed demonstrate that the magnet ic field strength B of FSRQs is different 
from and larger than that of BL 

Laci™, 

also consistent with our results. So the 
blazar sequence may be due to the different 7b and B of these sources. 



4.3. Implications for Blazar Sequence and Anti-Sequence 

[qi 

As shown in Fig. 1 of this paper and Fig. 1 of Ref.-^ 1 , the characteristics of spectral 
evolution at different states for a given source are opposite to the blazar sequence. 
In order to investigate the physical connections of the two phenomena, firstly, we 
need to define another physical parameter (Ry), the ratio of physical parameters 
in the high state to that in the low state for the twelve sources, where Y stands for 
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v s , Lboi, Lic/L syn , 7b, and U . Comparing the spectral evolution for a given source 
with the blazar sequence, we find the following: 

(i) According to the blazar sequence, the peak frequency of synchrotron radiation 
increases with decreasing bolometric luminosity. The ratio of bolometric lumi- 
nosity (i?L bol ) as a function of the ratio of peak frequency (R Us ) for the two 
states is presented in Fig. 6(a). It is found that for most of the sources in our 
sample the peak frequencies of synchrotron radiation move to the higher energy 
band in the high state than in the low state; at the same time the bolometric 
luminosities increase. 

(ii) According to the blazar sequence, the peak frequency of synchrotron radia- 
tion increases with decreasing ratio of the luminosities for the IC and syn- 
chrotron components. However, for most of the sources in our sample the values 
of L\c/L syn for high state become larger than that for low state, accompanied 
with the increase of synchrotron radiation peak frequency as shown in Fig. 6(b). 

(iii) A more "theoretical" characteristic for the blazar sequence is the anti- 
correlation between 7b and U in the comoving frame. However, as shown in 
Fig. 6(c), both 7b and U of the source in the high state become larger than 
that in the low state. 

As described above, both the blazar sequence from FSRQs to BL Lacs and anti- 
sequence of the individual sources from low state to high state are accompanied 
by an increase of 7b and a decrease of B, as shown in Fig. 3(a) and Fig. 5(c). 
Ghisellini & Tavecchio (2010) reported that the blazar sequence is linked to the 
different BH masses and accretion rates of different sourcesSSl. We propose here 
that the different states of the individual objects are also linked to the variations 
of the accretion rate for each source. So both the blazar sequence (change of source 
type) and the blazar anti-sequence of spectral variability of the individual objects 
(change of source state) are linked to the change of accretion rate. The flow chart 
illustrating how the accretion rate drives the blazar sequence and the anti-sequence 
is shown in Fig. 7, which is explained as follows: 

• Assuming the mass accretion rate M decreases, we can explain the blazar sequence 
(the left branch, (a-d) in the following) from FSRQs to BL Lacs and the blazar 
anti-sequence (the right branch, (e-h) in the following) for the individual sources 
from the low to high states. 

• (a) — Assuming equipartition between the magnetic energy and other forms of 
energies in the accretion disk and that the magnetic fields in the jets of blazars 
are carried over from the disk, then we expect B decreases. If 7b is determined 
by radiative cooling, then we expect 7b increases. 

• (b) — Statistically AGNs with smaller BH masses have lower M (in physical units) 
and thus lower bolometric luminosity Lboi- 

• (c) — The BLR exists only above a critical value of the accretion rate^ and its 
luminosity is well correlated with Lboi, though with a delay. It is well known 
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that the emission lines of BL Lacs are very weak, thus the second bump of their 
SEDs arc not contributed by IC/BLR (Leg)- Therefore, the values of L\c/L syn 
for FSRQs are larger that of BL Lacs as presented in Fig. 4(b). 

• (d) — From FSRQs to BL Lacs, J7 BLR and Ub (see Fig. 45(c)) decrease with 
the increasing 7b, resulting in an anti-correlation between U and 7b shown in 
Fig. 4(c). 

• (e) — For an individual source from the low to high state (M decreases), the 
electron energy -f c should increase with 7b (see Fig. 2(b) for larger 7b in high 
states) and thus the bolometric luminosity increases as shown in Fig. 6(a). This 
explains (i) above. 

• (f) — Ub decreases with M (as discussed above and see Fig. 3(a)), but U syn 
increases with 7b (for increasing L syn ). 

• (g) — Rli C /l s n increases for the BL Lacs in the high state as shown in Fig. 6(b), 
because Lie is "boosted" from L syn by 7b. For the two FSRQs, 6 is larger and B 
is smaller in high state than that in the low state, corresponding to larger U BLR 
and smaller Ub, so the ratio, Rh ic /L syn i is also larger than unity, as shown in 
Fig. 6(b). This explains (ii) above. 

• (h) — Less change oiUs and more increase of U syn result in R v i larger than unity 
for most of the BL Lacs, as shown in Fig. 6(c). This explains (Hi) above. 



5. Summary 

We have compiled the broadband SEDs of twelve TeV blazars that were simulta- 
neously or quasi-simultaneously observed with Fermi/LAT and other instruments 
from literature. Each of those sources has two broadband SEDs available, which are 
identified as a high or a low state according to its flux density at GeV/TeV band. 
We found that the syn+SSC model can well represent the observed SEDs for BL 
Lacs, whereas the EC/BLR contribution needs to be considered for explaining the 
observed SEDs of the two FSRQs. The magnetic field strength B of the two FSRQs 
is larger, but their Doppler factor 8 is smaller than that of the ten BL Lacs. Signifi- 
cant spectral shift to high energies accompanying with the flux increase is observed 
for each individual source, which seems opposite to the observational phenomenon 
of the blazar sequence. We refer this phenomenon to the blazar anti-sequence. How- 
ever, it is found that both the blazar sequence from FSRQs to BL Lacs and the 
anti-sequence of the individual sources from low to high states are accompanied by 
an increase of 7b and a decrease of B. We propose here that the blazar sequence 
and the anti-sequence of spectral variability of the individual sources are driven by 
decreasing accretion rate. A simple flow chart, which describes qualitatively how the 
accretion rate contributes to the blazar sequence and anti-sequence, is also given in 
this paper. 
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Fig. 1. The observed SEDs with model fitting for 3C 279 and PKS 1510-089. The data of high 
and low states are marked with blue and red symbols, respectively. The fitting parameters of 3C 
279 are p x = 2.1, p 2 = 4.4, 7 b = 294, At = 24 h, B = 4. G, <5 = 12 for the low state SED and 
pi = 2.2, p 2 = 4.46, 7 b = 526, At = 24 h, B = 2.65 G, 8 = 17 for the high state SED. The fitting 
parameters of PKS 1510-089 are pi = 1.1, p 2 = 3.8, 7 b = 631, At = 12 h, B = 2.1 G, 6 = 8.8 for 
the low state SED and pi = 1.9, p 2 = 3.2, 7 b = 400, At = 24 h, B = 1.15 G, 8 = 16 for the high 
state SED. 




Fig. 2. Panel a — Comparison of the peak frequencies u s (circles) and v c (triangles) between the 
high and low states. The solid line is the equality line. Panel b — Ratio R Jb as a function of the 
ratio R x TcV . The line is the best fitting line logR 7b = (-0.1 ± 0.3) + (0.49 ± 0.27) logRj TeV . 
The red symbols are the data for the FSRQs in our sample. 
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Fig. 3. The ratios of Ub, U h , L syn , and Lie as a function of the ratio of 7^ for the high state 
to the low state. The red symbols are the data for the FSRQs in our sample. 





Fig. 4. Panel a and Panel b — Bolometric luminosity Lbol and the ratio L c /L s as a function of v s . 
The data of high and low states are marked with blue and red symbols, respectively, and the stars 
are the data for two FSRQs. The best fit lines in Panel b are log(L c /L s ) = (4.0 ± 0.8) - (0.28 ± 
0.05) logi/ s for the low state data {dashed line) and log(L c /L s ) = (4.13±0.76) - (0.24±0.05) log u a 
for the high state data (solid line; excluding the source PKS 2005-489 in the high state). Panel 
c — 7b as a function of energy density U . The gray circles in Panel c are the sample data from 
Ghisellini et al. (2010). 
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Fig. 5. Panel a — The luminosity of BLRs (LBLa)as a function of synchrotron radiation peak 
frequency v s . The data are from Chen et al. (2009) . The back and the gray symbols indicate thermal- 
dominated and non-thermal-dominatcd FSRQs, respectively. Panel b, c — 7^ as the functions of 
r 2 and Ug. The data marked as circles are from Celotti & Ghisellini (2008). In the Panel c, black 
circles and light gray circles indicate BL Lacs and FSRQs, respectively, and the other symbols are 
the same as Fig. 4. 
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Fig. 6. Panel a, b — Ratios Rh hol and Rl c /l, as a function of the ratio R Vb for high and low 
states. Panel c — Ratio i? 7b as a function of the ratio R u i . The symbols are the same as Fig. 2. 
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Fig. 7. Flow chart on the mechanism behind the blazar sequence for different types of sources 
and anti-sequence for spectral variability of individual sources, both of which are driven by mass 
accretion rate M. The downward or upward arrows besides a parameter indicates decrease or 
increase of the parameter. Other arrows in the flow chart describe the underlying causal relations. 
See the text for detailed descriptions. 
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